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John Gemmef] 

Division of Applied Mathematics, Brown University, Providence, RI 02906, USA 

Eran SharorQ 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904 Israel 
Shankar VenkataramanfO 

Department of Mathematics, University of Arizona, Tucson, AZ 85721, USA 
(Dated: July 10, 2015) 

The edge of tom elastic sheets and growing leaves often form a hierarchical buckling pattern. Within non- 
Euclidean plate theory this complex morphology can be understood as low bending energy isometric immersions 
of hyperbolic Riemannian metrics. With this motivation we study the isometric immersion problem in strip and 
disk geometries. By finding explicit piecewise smooth solutions of hyperbolic Monge-Ampere equations on a 
strip, we show there exist periodic isometric immersions of hyperbolic surfaces in the small slope regime. We 
extend these solutions to exact isometric immersions through resummation of a formal asymptotic expansion. 
In the disc geometry we construct self-similar fractal-like isometric immersions for disks with constant negative 
curvature. The solutions in both the strip and disc geometry qualitatively resemble the patterns observed exper¬ 
imentally and numerically in torn elastic sheets, leaves and swelling hydrogels. For hyperbolic non-Euclidean 
sheets, complex wrinkling patterns are thus possible within the class of finite bending energy isometric immer¬ 
sions. Further, our results identify the key role of the degree of differentiability (regularity) of the isometric 
immersion in determining the global structure of a non-Euclidean elastic sheet in 3-space. 

PACS numbers: 02.40.-k, 68.90.+g, 87.10.Pq, 89.75.Kd 

The rippling patterns observed in torn elastic sheets (]J|3), 
leaves EH3 and swelling hydrogels 18-10] provide striking 
examples of periodic and self-similar patterns; see Fig. [T] 

Within the formalism of finite elasticity, it is understood that 
such patterns arise from the sheet buckling to relieve growth 
induced residual strains CEUH2. While numerical experi¬ 
ments set within this framework have been able to qualita¬ 
tively replicate these patterns mmm, there is no complete 
theoretical understanding of the mechanism behind the self¬ 
similar patterns. On the one hand, such complex, self simi¬ 
lar patterns are usually associated with “strongly frustrated” 
systems, e.g. elastic sheets with boundary conditions that 
preclude the possibility of relieving in plane strains asm 
or at alloy-alloy interfaces between distinct crystal structures 
ifTTl 181. On the other hand, many growth patterns generate 
residual in-plane strains which can be entirely relieved by the 
sheet forming part of a surface of revolution Q or a helix [16|. 

Given that generically this system is not strongly frustrated, 
why then do we observe self-similar buckling patterns? In 
this Letter we address this puzzle by showing that for a large 
class of growth profiles there exist periodic and self-similar 
deformations of the sheet with vanishing in-plane strain. The 
construction of these surfaces consists of gluing together local 
solutions of an isometric immersion problem along “lines of 
inflection” and “branch points” in such a manner that the re¬ 
sulting surface has finite bending energy. We propose that the 
sheet introduces these defects to locally reduce large bending 
content which necessarily results from the extrinsic geometry 
of the isometric immersions. 

One model of swelling thin elastic sheets is the non- 
Euclidean formalism of elasticity. This model posits that 



FIG. 1. Examples of periodic and self-similar wrinkling patterns 
in swelling thin elastic sheets, (a) Hydrogel disk with non-uniform 
swelling pattern, (b) Ornamental echeveria house plant, (c) Edge of 
a torn trash bag. 


growth permanently deforms the intrinsic distance between 
material coordinates which is encoded in a Riemannian metric 
g. By Gauss’s Theorema Egregium, g generates an intrinsic 
definition of Gaussian curvature K throughout the sheet fT9l 
in which points where K < 0 (K > 0) correspond to local 
swelling or growth (shrinkage or atrophy) l20l . The realized 
conformation of the sheet in three dimensional space is then 
one that minimizes bending energy over all isometric immer¬ 
sions of g EQ. 

For K uniformly negative, however, the extrinsic geome¬ 
try of the system imposes that with increasing domain size 
smooth isometries will develop singularities where one of the 
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principal curvatures diverges 1221 and for K = -1 these sin¬ 
gularities form curves - “singular edges” - [23 ]; see Fig. [2] 
Across the singular edge the bending energy diverges and 
hence with increasing domain size a different smooth isometry 
will be selected by energy minimization. The energetic price 
paid for switching to another isometry is a global increase of 
the bending content throughout the sheet; indeed for K = -1 
the maximum principal curvature of smooth isometries grows 
exponentially with domain size (24|. For the case of K = -1 it 
was shown that to mediate the exponentially growing bending 
content an energetically favorable alternative is for the sheet 
to adopt the configuration of an rc-fold periodic “monkey sad¬ 
dle” isometry ED These monkey saddles are not smooth 
across lines of inflection - every smooth isometry locally has 
the shape of a saddle - but nevertheless have finite bending 
energy. However, for a sheet with small but non-zero thick¬ 
ness, it has been shown from force and moment that these 
defects are in fact smoothed out by boundary layers yielding a 
smooth local minimizer of the elastic energy [25 j. That is, it is 
only in the vanishing thickness limit that these local minimiz- 
ers become non-smooth. The possibly non-smooth isometries 
selected by the energy in the vanishing thickness limit have 
recently been coined “asymptotic isometries” f26l . 



FIG. 2. Singular rim bounding surfaces of revolution with constant 
negative Gaussian curvature, e.g the tractricoid (a.k.a “the” pseudo¬ 
sphere). Near the singular edge one of the principal curvatures di¬ 
verges and the bending energy becomes infinite across the curve. 


To better understand the role the existence of the singular 
edge plays in the selection of possibly non-smooth isometries 
we now consider a strip geometry -oo<x<oo,0<y<oo 
with metric 


g = (l + e 2 f(y))dx 2 +dy 2 , (1) 

where / is a function satisfying f(y) > 0, f(y) —> 0 as y —> oo, 
f"(y) > 0 and 6 > 0. This metric corresponds to y depen¬ 
dent growth in the x direction localized near the y = 0 edge 
of the sheet and is a generalization of the metrics considered 
in EH30HSII2ZI. For £ 1 approximate isometries can be 

constructed by considering a Foppl - von Karman ansatz with 
order e 2 in-plane deformations and order e out-of-plane de¬ 
formation ew°(x,y). To order e 2 the solvability condition for 
an isometry is the following small slopes version of Gauss’s 


Theorema Egregium: 

9 n f”(y ) 

det (DV(v)) = -y, (2) 

where D 2 (w°(x,y)) denotes the Hessian of w°(x,y) m- 
The metric given by Eq. [T] admits exact smooth isome¬ 
tries in the form of surfaces of revolution □ or helices m, 
however as in the K = -1 case they have large bending con¬ 
tent throughout the domain. We now seek the analogue of 
the periodic monkey saddle solutions for this domain. We 
find that product solutions to Eq. ([2]) in the form w°(x,y) = 
2“ 2 \jj(x)(p(y) exist if for 0 < S < 1 we take 


[28l(\+6){\ +y/l)- 26Kl - S) 
jexp(-2y//) 


0<d< 1 
8 = 1 


(3) 


where / > 0 is the length scale of the swelling and \j/ satisfies 


i //' 2 ± k 25 \^\ 25 = i , 


(4) 


where k is a constant of integration. We can identify solu¬ 
tions if/(x) with the motion of a unit mass particle in the poten¬ 
tial V(i/s) = ±^(k\ifs\) 26 and with the positive branch we get 
periodic solutions with amplitude l/k and half wavelength 
Ak ~ k~ l . Note, however, for d ^ 1 ^ is not smooth across 
the lines of inflection if/ = 0; see Fig. [3] Nevertheless, for 
1 /4 < 8 < 1 it follows that if/ is square integrable and thus 
these deformations have finite bending energy. 



FIG. 3. Piecewise smooth small slope isometric immersion with out- 
of-plane displacement w°(x,y) = \f/{x)(p(y) for 8 = 1/2. The out- 
of-plane displacement vanishes along the green lines i.e. if/(x) = 0, 
and the periodic surface is created by odd periodically reflecting the 
smooth surface across these lines. 


For S = 1 we now explore the possibility of “lifting” 
the small slope isometries to full isometric immersions us¬ 
ing the formal asymptotic expansion in e. By solving the 
equations for an isometry order by order it follows that the 
order e 2a+l corrections to the out-of-plane displacement w = 
2 m) £ 2a+1 w a are of the form 

a 

w a (x,y) = k~ l e -( 2a+v >yl l ^ a n a ((k/) -1 ) cos((2^ + 1 )kx), 

n =0 

where a n j are polynomials. Now, the full series representa¬ 
tion w is only asymptotic in 6, i.e. it could diverge. Indeed, 
for increasing values of e - equivalently Gaussian curvature 
- we expect the surface to develop a singular edge, with the 
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first singularity appearing at (0,0), i.e. where |<9^w| attains 
its maximum value. If d 2 y w has a singularity of the form 
dyyW( 0, y ) ~ A(y + £ - £o) _( ^ +1) near (0,0) we can approximate 
£o and y 3 by the poles and residues of the Pade approximants 
to the logarithmic derivative d 2 y w( 0, 0)/d y w(0 ,0) l29l . In Fig. 
[4] we plot the Dlog Pade approximants for the case ( kl ) = 1 
and find that £o ~ 2.86 and = .500158. This strongly indi¬ 
cates that the principal curvature K y scales like s - 2 , where s 
is the arclength measured from the singular edge; a result that 
is consistent with the singular edge of the pseudosphere and 
other known hyperbolic surfaces of revolution. 



FIG. 4. [N, N + 1] Pade approximants to the logarthimic deriviative 
dy y w/dyW evaluated at (0,0). The inset figure is a plot of the aproxi- 
mate value of e 0 determined from the [N,N + 1] Pade approximants. 


The principal curvature K y diverges with increasing e for 
all values of k and /. This has consequences for modeling 
non-Euclidean elastic sheets. In particular minimizers of the 
small slope approximation to the elastic energy do not ade¬ 
quately approximate the full elastic energy. E.g, with 6 = 1 
and kl - 1 the solution w°(x,y) = k~ l e~ y/l cos(kx) is har¬ 
monic and hence a global minimizer of the bending energy 
over small slope isometries. However, in Fig. [5ja ) and 
Fig. [5jb) we present contour plots of the [9,9] Pade ap¬ 
proximants to Aw evaluated at (0,0) and the following proxy 
for the bending energy averaged over one half wavelength: 


b) indicates that for 


B[Aw] = A-' Xij /2 Jo° (Aw) 2 dydx. Fig. 
moderate values of 6, i.e. 6^2, the full elastic energy selects 
shorter wavelengths than predicted by the small slope theory, 
i.e. £ <$c 1. In particular, due to the existence of the singu¬ 
lar edge, it follows that for large values of e the wavelength 
must satisfy kl > 1 in order for the isometric immersion to 
exist. In particular, in a similar manner to the case K = -1 
to “go past” the singular edge it is necessary for this class of 
periodic isometries to refine their wavelength. The insets in 
Fig. [5ja) also indicate that the series solutions have a unique 
global wavelength and are smooth up to the singular edge. 

We now consider how periodic isometric immersions can 
introduce further defects in order to lower their bending con¬ 
tent. For simplicity we consider solutions to the Monge- 
Ampere equation det(Z) 2 w°) = -1, the small slopes approx¬ 



FIG. 5. (a) Contour plot of the [9,9] Pade approximants to Aw evalu¬ 
ated at (0,0) as a function of ( kl)~ l and e for the specific value 6=1. 
The dashed curve indicates values at which the singular curve first 
appears at (0,0). (b) Contour plot of the bending energy B [Aw] as a 
function of ( kl)~ l and e for the specific value 6=1 using the [9,9] 
Pade approximants to Aw. 


imation to K = -1, on the disk x 2 + y 2 < R 2 . In ref¬ 
erence ci, the n-fold periodic solutions to this equation 
were formed by taking odd periodic reflections of the solu¬ 
tion w„(x,y) = y(x - cot(7r /n)y) about the line of inflection 
y = tan(7r/ft)x; see Fig. [6ja-b). We adopt the terminology in 
(23l and refer to points where lines of inflection intersect - 
in this case x = y = 0 - as branch points. They correspond 
to branch points of the map (v,y) i-> Vw, or equivalently to 
branch points for the Gauss normal map [ 30]; see Fig. [6jc-d). 

The key to the above construction is that any quadratic 
surface is ruled by two sets of parallel asymptotic lines and 
































4 


(a) (b) 





FIG. 6. (a-b) Small slope isometric immersions wj}(x, y) and w®(x, y) 
for constant Gaussian curvature K = - 1. j) is constructed by 
taking odd periodic reflections of the piece of w®(x,y) bounded be¬ 
tween the green lines. The mesh on both of these surfaces correspond 
to their asymptotic lines, (c-d) Direction of the gradient Vw along 
circles centered at the origin. The regular saddle in (a) corresponds 
to a gradient field with winding number -1, so the gradient map is 1 
to 1. The 4-saddle in (b) has winding number -3, so the gradient map 
is a 3 sheeted covering near the origin. 



FIG. 7. Finite bending energy solutions to the Monge-Ampere equa¬ 
tion [w°,w°] = -1. (a) Three sub wrinkle solution created by in¬ 
serting three rotated and translated copies of the solution Wg(v,y) = 
y(x- V3 y) onto the solution \^(x,y) = xy at a branch point, (b) Nine 
subwrinkle solution created by inserting nine copies of w° l2 (x,y) = 
y(x - (2 + V3)y) at three branch points added onto the three sub¬ 
wrinkle solution, (c) Extension of the nine subwrinkle solution to 
the full circular domain, (d) The nine subwrinkle solution mapped to 
the strip geometry by a conformal map. 


to retain continuity of the tangent plane the lines of inflec¬ 
tion are introduced along asymptotic lines; see Fig. 0b). 
Replicating the above process at other points, multiple branch 
points can be introduced on the surface not just at the ori¬ 
gin. For example, for w Q 2 (x,y) = xy a branch point b\± can 
be added by removing the sector x,y > 1 / V2 and in this re¬ 
gion fitting three appropriately rotated and translated copies of 
w%x,y) = y(x- V3j); see Figma). Three more branch points 
Z? 2 ,i? ^ 2 , 2 , ^ 2,3 can be added along rays emanating from b\^\ 
that bisect the lines of inflection; see Fig^b). The solution 
can be extended by odd periodic reflections to give a small- 
slopes isometric immersion of the unit disk with K = -1. To 
illustrate the wrinkling behavior near the edge we map w to a 
strip geometry through a conformal map h[x + iy ] = w[e x+iy ]; 
see Figs. 0c-d). 

The existence of isometric immersions with branch points 
also has implications to the modeling of non-Euclidean elastic 
sheets. As for the strip with S = 1, the solution w^(x,y) is 
harmonic yet the extension of w^{x,y) to an exact isometric 
immersion has divergent bending energy for R > 1.25 with 
the bending content concentrated near the singular point x = 
y « 1.25/ V2 [241. We can isometrically immerse disks with 
larger R by a global refinement of the wavelength i.e taking 
n > 2. These solutions increase the bending energy globally. 
An energetically favorable alternative might be to introduce 
a branch point in the n = 2 solution near the singular point, 
and locally refine the wavelength instead. Indeed, numerics 
for S = 1/3 in the strip geometry indicate that, even within the 


small slopes approximation, localized self similar wrinkling 
profiles may be energetically preferred over global refinement 
of the wavelength GJQJ). 

The existence of finite bending energy isometric immer¬ 
sions for hyperbolic free sheets ensures that the per unit thick¬ 
ness energy scales like t 2 . This is in contrast to crumpled 
sheets which have an energy scale t 5/3 which is intermedi¬ 
ate between the stretching and bending energies [31, 32]. 
Furthermore, our results show that there are multiple low- 
energy states as finite bending energy isometric immersions 
can be constructed by appropriately gluing together low en¬ 
ergy building blocks using lines of inflection and branch 
points in a variety of ways. The energy barriers between dis¬ 
tinct low energy states are small, ~ t 2 , and these sheets are 
thus “floppy”. Consequently, thin elastic sheets are easily de¬ 
formed by weak stresses, and the pattern selected may be sen¬ 
sitive to the dynamics of the swelling process, experimental 
imperfections, or other external forces. A key first step for 
better understanding the buckling patterns of hyperbolic free 
sheets would be to analyze how the optimal elastic energy of 
the two types of singularities - lines of inflection and branch 
points - scales with the various length scales in the problem: 
k ~ l , l and t. 
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